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Pygmy and Giant Dipole Resonance (FDR and GDR) in Ni isotopes have been investigated by 
Coulomb excitation in the framework of the Isospin-dependent Quantum Molecular Dynamics model 
(IQMD). The spectra of 7 rays are calculated and the peak energy, the strength and Full Width 
at Half Maximum (FWHM) of GDR and PDR have been extracted. Their sensitivities to nuclear 
equation of state, especially to its symmetry energy term are also explored. By a comparison with 
the other mean-field calculations, we obtain the reasonable values for symmetry energy and its slope 
parameter at saturation, which gives an important constrain for IQMD model. In addition, we also 
studied the neutron excess dependence of GDR and PDR parameters for Ni isotopes and found 
that the energy- weighted sum rule (EWSR) PDRmi / GDRmi % increases linearly with the neutron 
excess. 

PACS numbers: 25.75.Dw, 25.70.Do, 25.70.Ef 
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I. INTRODUCTION 



Intermediate energy heavy-ion colhsons provide a 
chance to study reaction machanism and the propertis 
of hot nuclei, such as fragmentation and hquid gas phase 
transition etc In a relative lovi^ excitation region 

of nuclei, dipole resonance has attracted much attention 
experimentally and theoretically in the past few decades 
H • Two kind of resonances have been evidenced, one 
is Pygmy Dipole Resonance (PDR) and another is Giant 
Dipole Resonance (GDR). In contrast to GDR, which can 
be considered as the oscillation between non-deformed, 
incompressible proton and neutron spheres, PDR can be 
considered as the oscillation between the weakly-bound 
neutron skin and the isospin neutral proton- neutron core. 
From a theoretical point of view, the presence of this 
low-lying strength PDR is predicted by almost all micro- 
scopic models, ranging from Hartree-Fock plus random 
phase approximation (RPA) with Skyrme interactions to 
relativistic Hartree-Bogoliubov plus relativistic quasipar- 
ticle RPA 0. Experiments of PDR jj-fl^ in neutron- 
rich nuclei and microscopic models which have been ap- 
plied to investigate PDR [l5l - [2l| have shown that the 
PDR could have a pronounced relationship with neutron- 
capture rates in the 7-process, nucleosynthesis, the radia- 
tive neutron-capture cross section on neutron-rich nuclei, 
and the photodisintegration of ultra-high energy cosmic 
rays. The properties of the PDR in stable nuclei have 
been studied extensively for different neutron and pro- 
ton shell closures with the (7,7') reactions. However, 
studies for neutron-rich isotopes, where the PDR should 
be enhanced, were still less, especially in the framework 
of molecular-dynamics type model. Recently, the influ- 



*Author to whom all correspondence should be addressed. Email: 
ygma@sinap.ac.cn 



ence of symmetry energy on PDR has been also checked 
p^ ? ]. For example, in Ref. [? ], they studied isovcctor 
giant and pygmy dipole resonances in even-even Ni iso- 
topes within the framework of a fully consistent relativis- 
tic random-phase approximation built on the relativistic 
mean field ground state and found that the centroid en- 
ergy of the isovector pygmy resonance is insensitive to 
the density dependence of the symmetry energy. Ref. (22j 
used a semiclassical Landau- Vlasov approach to investi- 
gate the energy centroid associated with the PDR and 
found that PDR is insensitive to the symmetry energy 
term of EOS. Considering symmetry energy is one of in- 
teresting topics in current heavy-ion physics community, 
we would also like to check this sensitivity in the present 
work. 



In this work, we apply the Isospin-dependent Quantum 
Molecular Dynamics model (IQMD) to investigate the 
properties of PDR. In our previous work, the same model 
has been successfully applied into the GDR calculation. 
Here we plan to extend the model to calculate PDR in 
the same framework. The molecular dynamics model is 
a kind of Monte-Carlo transport model, which has been 
extensively applied in heavy-ion collision dynamics. In 
the model, the physical parameters can be controlled and 
then their respective effects on PDR and GDR can be 
addressed. 



The paper is organized as follows. Sec. II gives a brief 
introduction of IQMD model as well as the formalism 
to calculation GDR and PDR in the IQMD framework. 
Results and discussions are presented in Sec. Ill where 
different parameter dependencies of GDR and PDR are 
shown. In particular, the equation of state, symmetry 
energy and isospin effects are presented and discussed. 
Finally the summary is given in Sec. IV. 
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II. MODEL AND FORMALISM 

The isospin-dependent quantum molecular dynamics 
(IQMD) model is based on QMD model [H. In IQMD 
model, the mean field is given by: U{p) = U^^^ + U'~^°'"-'- + 

jjYuk _|_ jjsym _|_ jjMDI ^ where U^^y ^ jjCoul ^ jjYuk ^ Jjsym 

and JJ^^^^ is the Skyrme potential, Coulomb poten- 
tial, Yukawa potentia, symmetry potential interaction, 
and the momentum dependent interaction (MDI) , respec- 
tively. The Skyrme potential can be presented as follows: 



(1) 



where po = 0.16/ fm^ (the saturation nuclear density) 
and p is the nuclear density. In the equation, different 
[a, P, 7] represents different kinds of equations of state 
(EOS). The parameters a, /3 and 7 are given in Table H] 
In the Table, S(M) represents the soft EOS (with MDI), 
with an incompressibility of K~2QQ MeV, while II(M) 
for the hard EOS (with MDI), with an incompressibility 
of A'=380 MeV. 



TABLE I: The parameters a, P and 7 for the different EOS 



EOS 


K 


a 


P 


7 




(MeV) 


(MeV) 


(MeV) 


(MeV) 


S 


200 


-356 


303 


7/6 


SM 


200 


-390.1 


320.3 


1.14 


H 


380 


-124 


70.5 


2 


HM 


380 


-129.2 


59.4 


2.09 



jjCoui^ jjYuk^ jjsyni jjMDi ^^^^ exprcsscd as 
follows, respectively: 
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Coul 



e 

T 



E 



1 



(47rL)(3/2) 



exp[ 



4L 



(2) 



Here, 



r.; - r,- 



represents the relative distance of 



two nucleons, and the L is the so-called Gaussian wave- 
packet width for nucleons. For simplicity, here we choose 
a constant wave-packet width, i.e. L= 2.16 fm^. We 
noted that there was some discussions on the affect of 
the width on the dynamical results, eg. flow, multifrag- 
mentation, pion and kaon production etc [23 |. However, 
we are treating the Coulomb excitation of the projectile 
in which the dynamical effect is relatively smaller. 



X [exp{mrij)erfc{\/Lm — i^ij/y/AL) (3) 
-exp{mrtj)erfc{\/Lm + r.j/V^I)], 

where Fy=0.0024 GcV and m = 0.83. 



TTSym _ C^sym \ ^ ^ 

- 2po ^"-"^■^(4.L)3/2 



exp[-^^^^], (4) 



where Cgym is the symmetry energy coefficient, is the 
zth component of the isospin degree of freedom for the 



nucleon, which equals 1 or -1 for neutron or proton, re- 
spectively. 



U 



MDI 



Po 



Po 



(5) 



Here, S =1.57 MeV and e =500 c^GeV^. More details 
of the QMD model can be found in [ll, IH-il . 

The dipolc moment of GDR in coordinator space 
[DRcDRit)] and momentum space [DKcDRit)] is respec- 
tively defined as follows [H, [s^ : 



DRoDRit) 
DKGDRit) 



NZ 



NZ ^ Pz{t) PN{t) ^ 
Ah 



z 



N 



(6) 
(7) 



where Rz (t) and Rn (t) are the center of mass of protons 
and neutrons in coordinator space, respectively; Pz{i) 
and Pn (t) are the center of mass of protons and neutrons 
in momentum space, respectively. 

By the Fourier transformation of the second derivative 
of DRgdrX^) respect to time, i.e. 



DR"{uo)= / DR"{t)e"^'dt, 

'to 



(8) 



we can get the the 7 ray emission probability for energy 



Ej = fiuj as follow: 



dP 



2e^ 



\DR"[ur 



(9) 



Similarly, we can write the dipole moment of FDR in 
coordinator space and momentum space as follows: 



DRpDRit) = -^[Rc{t) - RNAt)] 



DK 



PDR 



A 

N.,Z Pcjt) 
Ah ^ 2Z 



(10) 

(11) 



where Rc{t) and RN^{t) are the center of mass of 
isospin-symmetric core and valence neutrons in coordi- 
nator space, respectively; Pz{t) and Pwit) are the center 
of mass of isospin-symmetric core and valence neutrons 
in momentum space, respectively. Then we can also get 
the 7 ray emission probability of FDR. Here, we choose 
neutrons with the farthest distances from the CM of all 
nucleons as the valence neutrons in the initial state. 

The fraction of the energy-weighted sum rule (EWSR) 
contained in the FDR relative to that located in the GDR 
region can be written as follows 



^ ^^5^ X 100%, 



(12) 



GDR„i^ 

where GDR^^ and PDR,n, is EWSR of GDR and PDR, 



3 



respectively: 



,p 

GDR„,, - Y,i — )GDRAE X E, 



E'r, 



dp 



(13) 
(14) 



III. RESULTS AND DISCUSSIONS 

Now we move on the calculations and discussions. As 
a first step of this work, we need to check the stability of 
the nuclei [l^]- Here we select ^*Ni+^^^Au as an example, 
where the Soft EOS with MDI has been used. Figure [1] 
displays the binding energies and root mean square radii 
of the cold projectile (^^Ni) and target (^^^Au) nuclei as 
a function of time till 400 fm/c. It is clearly seen that 
both variables keep very stable during the whole time 
evolution, which assures that we can proceed with the 
following calculations on the GDR and PDR. 




50 100 150 200 250 300 350 400 
t (fm/c) 



FIG. 1: (Color online) Stability check for the projectile C^^Ni) 
and target (^^^Au). (a) time evolution of the binding energies; 
(b) time evolution of rms radii. 

After the stability check, we shall compare the cal- 
culated results with the experimental data to demon- 
strate the reliability of our calculation. The calculation 
parameters are as follows: incident energy (i?™) is 600 



MeV/nucleon, impact parameter (6) is 24fm, the Soft 
EOS with MDI, and C^ym is 32MeV. The PDR resuh of 
^^Ni is shown in Figure [2] where the symbol with error 
bar is for the experimental data from Ref . [ll| , the line is 
our calculated result. We can see that the peak energy of 
our calculated result has good agreement with the data 
even though the FWHM of our calculated result is a little 
larger than the data. For GDR calculations, our previous 
IQMD calculations have shown a good agreement with 
the data (soj . After this comparison, we perform a sys- 
tematic calculation to explore the variable dependencies 
of GDR and PDR. 




9 10 11 12 13 14 15 16 17 18 
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FIG. 2: (Color online) The calculated PDR result of ''^Ni 
compared with the experimental data. The blue circle with 
error bar is the experimental data in Ref. the black line 
is the calculated result. 

Firstly, the sensitivities of GDR and PDR of ^^Ni in the 
Coulomb excitation of ^^^Au to some parameters have 
been studied. From the calculated results, we can extract 
the peak energy, strength and FWHM of GDR or PDR 
by Gaussian fitting to the 7 emission spectrum. Figure [3] 
shows the time evolution of dipole moment for GDR and 
PDR at different incident energies. With an increasing of 
energy, both amplitude and frequency of GDR and PDR 
oscillations become smaller, this could be attributable to 
the weakening of the effective Coulomb field in which the 
projectile interacts with the target when the projectile 
has faster velocity. Actually this can be seen from the 
inverse of the distance (l/Rpx) between the CM of pro- 
tons of projectile and CM of protons of target. Figure |4] 
shows 1/RpT as a function of time in different incident 
energies. If we roughly treat the projectile and target as 
point charged particles, then the Coulomb interaction is 
proportional to the time integration of l/Rpx- There- 
fore, the Coulomb effect becomes smaller when the in- 
cident energy is larger and then the projectile gets less 
excited. 

On the other hand, in Figure [3] pygmy dipole oscillation 
shows a much smaller amplitude and a slightly smaller 
frequency in comparison with the GDR case. These lead 
to the behavior of the incident energy dependence of the 
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FIG. 3: (Color online) Time evolution of dipole moment for 
GDR (a) and PDR (b) at different incident energies. The 
meaning of lines are illustrated in the insert. In calculations, 
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MDI. 
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FIG. 5: (Color online) Incident energy dependence of dipole 
moments for ®*Ni. The upper panel for GDR and the bottom 
one for PDR. In calculations, we use b — 24 fm, Csym = 
32MeV, and the Soft EOS without MDI. 
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FIG. 4: (Color online) Time evolution of dipole moment for 
GDR and PDR at different incident energies. The meaning 
of lines are illustrated in the insert. In calculations, we use 6 



= 24 fm, C^ 



=32MeV, and the Soft EOS without MDI. 
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FIG. 6: (Color online) EOS dependencies of GDR (left pan- 
els) and PDR (right panels) parameters for ^*Ni. In calcula- 



tions, we use Ein- 
32MeV. 



600 MeV/nucleon, 6= 24 fm, and C^ 



GDR and PDR parameters as shown in Figure [5] . With 
the increase of incident energy, both the peak energies 
and strength for GDR and PDR drop. For the FWHM, it 
clearly shows an increasing behavior for GDR. Similarly, 
FWHM of PDR displays an increasing behavior versus 
energy within large errors. 

Figure H] shows EOS dependence of the dipole res- 
onance parameters. Compared to the Soft EOS with 
or without momentum dependent interaction, the calcu- 
lated results with the Hard EOS arc generally bigger. On 
the other hand, compared to the EOS without MDI, all 
the peak energy, strength and FWHM in the EOS with 
MDI are also larger. In this sense, MDI has similar effect 
to make the EOS harder. From this figure, we learn that 
diploc 7-cmission is very sensitive to the stiffness of EOS. 



Many works have demonstrated that the symmetry en- 
ergy plays an important role in understanding the mech- 
anisms of many exotic phenomena in nuclear physics and 
astrophysics. Dipole oscillations between all neutrons 
and all protons or between the excess neutrons and a 
core composed of an equal number of protons and neu- 
trons involve asymmetry of neutron and proton, which 
could be infiuenced by the symmetry energy term of the 
EOS. To this end, we change the value of Csym in IQMD 
model to address its effect on GDR and PDR. In order to 
see a clear trend of GDR and PDR toward the symmetry 
energy coefficient, we choose a larger range of Csym in 
the calculations, i.e. from 16 - 64 MeV. In between, a 
value of Csym around 36 MeV has been thought reason- 
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able to describe the property of ground state of nuclei in 
many previous studies. Figure [7] shows time evolution of 
dipole moments for GDR (upper panel) and PDR (lower 
panel) in different symmetry coefficients Csym- It shows 
that the frequency of GDR oscillations becomes faster 
and the amplitude tends to be smaller with the increas- 
ing of Csym, which leads to an increasing of peak energy 
and decreasing of strength as shown in next figure. On 
the contrary, the frequency of PDR oscillations becomes 
slightly slower, but the amplitude tends to be smaller as 
well as GDR with the increasing of Csym, which induces 
a decreasing of peak energy as well as strength as shown 
in next figure. 
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FIG. 7: (Color online)Cst;m dependence of of dipole moments 
for ^*Ni. The upper panel for GDR and the bottom one for 
PDR. In calculations, we use Ein= 600 MeV/nucleon, b — 24 
fm, and the Soft EOS without MDI. 

Figure [5] shows the calculated results of PDR and GDR 
parameters with the different Csym parameters. With 
the increases of Csym, the peak energies of GDR show 
a linear increase, but those of PDR shows a little de- 
crease. However, considering the obvious different scale 
in y-axis for GDR and PDR, peak energy of PDR has 
only about 4% decrease in contrast to 25% increase for 
GDR. For strength, both show a decreasing behavior and 
for FWHM, both display an increasing trend. In some 
previous studies, it was known that the restoring force 
of isovector GDR is proportional to the symmetry en- 
ergy of nuclear matter, which makes the peak energies 
of GDR pushed to higher energies with stronger symme- 
try energy. By contrast, PDR originates mainly from the 
vibrations of a few valence neutrons which against the 
isospin neutral core. Those valence neutrons are located 
in the very low density exterior region, where the sym- 
metry energy is not affected as much. In general, the 
main difference by C'sym is observed in the interior re- 
gion where is responsible for GDR, which has very little 
or weak infiuence on the subtle properties of PDR. Sim- 
ilar insensitivity of PDR to symmetry energy was also 
found in some literatures [22l ? ]. 

It will be very interesting if we can pin down some in- 




FIG. 8: (Color online) Csym dependence of GDR (left panels) 
and PDR (right panels) parameters for ^*Ni. Symbols present 
for calculations and curves are linear fits for guiding the eyes. 
In calculations, we use Ei„= 600 MeV/nucleon, b — 24 fm, 
and the Soft EOS without MDI. 
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FIG. 9: (Color online) 

{PDR^JiPDRm, + GDR,^,)%) 



Calculated EWSR 
for ^*Ni with differ- 



ent Esym in our IQMD model (open stars) together with 
a bunch of mean- field model calculations (solid circles), 
which is taken from Ref. [3l|]. Line is a fit to mean- field 
model calculations. In the figure, A, B, C and D represent 
mean-field calculations with SkMP, SkRs, SkGs, and SK255 
parametrization. For details, see Ref. [33|. In figure, the 
shaded area represents 15% uncertainty around the point A 
which is just on fitted line. 



formation on the symmetry energy, such as the Esymipo) 
and the derivative of the symmetry energy at saturation 
(Ls) by the dipole resonance calculation. In our IQMD 
model, the symmetry energy reads 



Esym = I2.5(-^)2/3 

Po 



Csym / P \ 

2 Po 



(15) 
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where the first term is the kinetic energy contribution and 
the second term is potential energy term. In this case, 
the slope parameter of the symmetry energy at satura- 
tion, Ls = 3Csy,„/2-|-25. With the various Csym, we can 
obtain a plot of EWSR P D R^J {P D + GDR„,,)% 
versus Ls as shown in Figure [9l To compare with other 
model calculations, different mean-filed calculation re- 
for ^^Ni are displayed in the figure. From the 



(a) GDR 



suits 31 



figure, we observe a weak anticorrelation of EWSR vs Ls 
in IQMD calculations. However, the IQMD results are 
quite similar to the mean-filed calculations with SkMP, 
SkRs, SkGs, and SK255 parameterizations, i.e. the corre- 
sponding A, B, C and D points [3l|. The anticorrelation 
is not so reasonable in comparison with the overall posi- 
tive correlation trend which has been constructed by the 
mean-field calculations from Ls ~ 10-120 MeV. This may 
originate from the unreal symmetry energy coefficient 
when Csym is selected too large or too smaller. Assum- 
ing the fitting line is correct for the relationship of EWSR 
versus Lg and allow 15% uncertainty, we can obtain that 
the reasonable Ls values from 60 to 81 MeV for IQMD 
calculations (see the shaded area), which leads to Cgym ~ 
23.3 - 37.3 MeV , and gives Esym{po) ~ 24.2 - 31.2 MeV. 
The above values for IQMD are in good a gree ment with 
other calculations. For instance, in Ref. [31[, they got 
the weighted average, Ls = 64.8 ± 15.7 MeV, and the 
deduced best value of Esym{po) = 32.3 ± 1.3MeV. In 
Ref. im they give the value of Esym{po) = 32.0 ± 1.8 
MeV and in Ref. [s^ the ranges obtained of Esym{po) is 
30.2 - 33.8McV. So far, the above Esym{po) and L pa- 
rameters give a useful constraint for the IQMD model. 

We also studied the systematic evolution of GDR and 
PDR parameters of even-even Ni isotopes from ^^Ni to 
^^Ni. Figure [TU] shows the distributions of isovector 
dipole strength for the Ni isotopes. In the figure, the 
GDR peak is around 12 MeV and that of PDR is around 
10 MeV. The PDR strengths increase as the neutron 
number increases. 

Figure [TT] shows the extracted GDR and PDR param- 
eters for Ni isotopes from Figure [TOl With the increases 
of mass number A or neutron excess, the peak energies of 
both GDR and PDR decrease. This behavior reflects that 
a larger neutron excess results in lower GDR and PDR 
excitation energy. In the region beyond A — 72, the slope 
of the PDR peak energies becomes steeper than those of 
lower mass nuclei, because the neutrons in outer orbital 
are more loosely bound and thus the restoring force in the 
oscillation of the skin against the core becomes weaker. 
The strength of GDR is shown to seemingly decrease de- 
creasing for larger mass nuclei but it increases for PDR, 
the FWHM of both GDR and PDR show increasing be- 
havior with A, except for smaller A for GDR. 

Finally, the fraction of the EWSR toi% is also ex- 
tracted. Figure [T^] shows the calculated results which 
clearly illustrate that with the increases of mass number 
A or neutron-skin thickness, PDR component increases 
almost linearly. Similar phenomenon has been also re- 
ported earlier for some isotopes with random phase ap- 
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FIG. 10: (Color online) Mass number dependence of GDR (a) 
and PDR spectra (b) for Ni isotopes. In calculations, we use 
Ei„= 100 MeV/nucleon, 6= 24fm, Csym =32 MeV, and the 
Soft EOS without MDI. 



proximation phenomenological approach [20|, |31|, |3J, [3 



IV. SUMMARY 

In summary, we have applied the IQMD model to 
study Giant and Pygmy Dipole Resonance in Ni isotopes 
by the Coulomb excitation. Similar to the method to 
calculate Giant Dipole Resonance in our previous IQMD 
calculation, we extend it to calculate PDR in the same 
model. We first showed very stable initial projectile 
and target can be obtained with the soft momcnteum- 
dependent EOS, which assures the performance of GDR 
and RDR calculations. After we got a satisfied calculated 
PDR result for ^^Ni which agrees with the experimental 
data very well, then we performed a systematic calcula- 
tion for GDR and PDR with different beam energy, EOS 
and symmetry energy parameters for ^^Ni-f-^^^Au situa- 
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FIG. 11: (Color online) Mass number dependence of Ni iso- 
topes of GDR (left panels) and PDR (right panels) param- 
eters. In calculations, we use Ein~ 100 MeV/nucleon, b— 
24fm, Csym =32 MeV, and the Soft EOS without MDI. 
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FIG. 12: (Color online) Calculated EWSR 

PDRm^/GDRm^% for Ni isotopes (solid circles) together 
with a linear fit (line). 



tion. It is found that the peak energy and strength for 
both GDR and PDR show a decreasing trend with the 
incident energy, which can be understood by the exci- 
tation extent due to Coulomb field in different energies. 
When the equation of state becomes harder or momen- 
tum dependent interaction is taken into account, the in- 
creasing behaviors of peak energy, strength and FWHM 
emerge for both GDR and PDR. Concerning the sym- 
metry term of EOS, it shows strong positive correlation 
with the GDR's peak energy and FWHM, but it plays a 
relative weak role on PDR. By fitting to other mean-field 
calculations of EWSR vs the slope parameter at satura- 
tion Ls and assuming 15% uncertainty of Lg, i.e. Lg = 
60 -81 MeV, we can obtain the range of symmetry energy 
coefficient Csym ~ 23.3 - 37.3 MeV and symmetry energy 
at saturation is around 24.2 - 31.2 MeV and in the IQMD 
model, which gives a constraint for the QMD model. Fi- 
nally, mass number or isospin dependence of GDR and 
PDR for the Ni isotopes is presented. The result shows 
that the peak energy of both GDR and PDR decreases 
with mass number, and the FWHM of them increases. 
For the strength, PDR shows an increasing behavior but 
GDR seemingly decrease with mass of isotopes. For the 
fraction of the EWSR in Ni isotopes, it shows a linear 
rise-up with the increases of mass number A or neutron 
excess. 
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